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Abstract. We consider multihadron production processes in different types of collisions in the 
framework of the picture based on dissipating energy of participants and their types. In particular, 
the similarities of such bulk observables like the charged particle mean multiplicity and the pseu- 
dorapidity density at midrapidity measured in nucleus-nucleus, (anti)pro ton-pro ton and electron- 
positron interactions are analysed. Within the description proposed a good agreement with the 
measurements in a wide range of nuclear collision energies from AGS to RHIC is obtained. The 
predictions up to the LHC energies are made and compared to experimental extrapolations. 
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1. High densities and temperatures of nuclear matter reached at RHIC provide us with 
an exceptional opportunity to investigate the matter at extreme conditions. Bulk observ- 
ables such as multiplicity and particle densities (spectra) being sensitive to the dynamics 
of strong interactions, are of fundamental interest. Recent measurements at RHIC re- 
vealed striking evidences in the hadron production process including similarity in such 
basic observables like the mean multiplicity and the midrapidity density measured in 
complex ultra-relativistic nucleus-nucleus (AA) collisions vs. those obtained in rela- 
tively "elementary" e^e^ interactions at the same centre-of-mass (cm.) energy when 
number of participants ("wounded" nucleons IH] in AA collisions) are taken into ac- 
count ||2|,|3||]. The observation is shown to be independent of the cm. energy per nucleon 
y/^NN = 19.6 GeV to 200 GeV. Assuming similar mechanisms of hadron production in 
both types of interactions which then depends only on the amount of energy transformed 
into particles produced, one would expect the same value of the observables to be ob- 
tained in hadron-hadron collisions at close can. energies. However, this is not the case: 
comparing measurements in hadronic data [4, 5] to the findings at RHIC, one obtains 
|2,^XI quite lower values in hadron-hadron collisions. In the meantime, the RHIC dAu 
data at a/^nn ^ 200 GeV unambiguously point to the values of the mean multiplicity 
from pp data [2]. Moreover, recent CuCu RHIC data show no changes in the values of 
the bulk variables compared to those from AuAu collisions when properly normalised 
to the number of participants [8, 91. 

The observations made earlier|0] and the recent ones |0] can be understood in the 
franework of a description proposed recently by us ifiol] and considered here. This 
description is based on a picture when the whole process of a collision is interpreted as 
the expansion and break-up into particles of an initial state, in which the whole available 



energy is assumed to be concentrated in a small Lorentz-contracted volume. There are 
no any restrictions due to the conservation of quantum numbers besides energy and 
momentum constraints allowing therefore to link the amount of energy deposited in 
the collision zone and features of bulk variables in different reactions. This description 
resembles the Landau hydrodynamical approach to multiparticle production [TT] which 
has been found to give good description of the mean multiplicity AA, pp, e+e , vp data 
|12, 13] as well as of pseudorapidity distributions at RHIC |@]. 

As soon as a collision of two Lorentz-contracted particles leads to the full thermaliza- 
tion of the system before extension, one can assume that the production of secondaries 
is defined by the fraction of participants energy deposited in the volume of the system 
at the collision moment. This implies that there is a difference between results of col- 
lisions of structureless particles like electron and composite particles like proton, the 
latter considered to be built of constituents. Indeed, in composite particle collisions not 
all the constituents deposit their energy when they form the Lorentz-contracted volume 
of the thermalized initial state. As a result, the leading particles [14], formed out of those 
constituents which are not trapped in the interaction volume, carry away a part of en- 
ergy. Meantime, colliding structureless particles are ultimately stopped as a whole in the 
initial state of the thermalized collision zone depositing their total energy in the Lorentz- 
contracted volume and this energy is wholly available for production of secondaries. 

We consider a single nucleon as a superposition of three constituent quarks due to 
the additive quark picture fis"]. In this picture, most often only one quark from each 
nucleon contributes to the interaction with other quarks being spectators. Thus, the initial 
thermalized state is pumped in only by the energy of the interacting single quark pair 
and, so, only 1/3 of the entire nucleon energy is available for production of secondaries. 
Therefore, one expects that the resulting bulk variables like the multiplicity and rapidity 
distributions should show identical features in pp collisions at the cm. energy ^/s^ 




and in e+e^ interactions at the cm. energy ~ Note that for the mean 

multiplicity, a similar behaviour was found in the beginning of LEP activity [16]. 

In AA collisions, more than one quark per nucleon interacts due to the large size of 
nucleus and the long path of interactions inside the nucleus. In central AA collisions, a 
contribution of constituent quarks rather than par ticip ating nucleons seem to determine 
the properties of produced particle distributions iTltIi . In headon collisions, the density 
of matter is almost saturated, so that all three constituent quarks from each nucleon may 
participate nearly simultaneously in collision depositing their energy coherently into the 
thermalized zone. Therefore, in the headon AA interactions at ^/¥mn the bulk variables 
are expected to have the values similar to those from pp collisions at y^s^ — 3 ^/^nn- 




This makes the most central collisions of nuclei akin to e+e collisions at y/s^ ~ ^/JMN 
in sense of the resulting bulk variables. 

2. According to our consideration, in Fig.[ll we compare the cm. energy dependence 
of the mean multiplicity in AA and e+e interactions to that in pp/pp collisions at 
\Ann = = from a few GeV to 200 GeV. For y/s^ > M^o, we give the 

multiplicities averaged |31] from the recent LEP data at -s/s^ =130 GeV and 200 GeV: 
23.35 ± 0.20 ± 0. 10 and 27.62 ± 0. 1 1 ± 0. 16. Figure shows also the mean multiplicity fit 
to pp/pp data [20] and the 3NL0 pQCD |33] ALEPH fit to e+e data L27]. 

From Fig. [Hone sees that the pp/pp data are very close to the e+e" data at y/I^ = 
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FIGURE 1. The charged particle mean muhipHcity Nch per participant pair as a function of the 

cm. energy. The soHd and combined symbols show the multiplicity values from: most central heavy-ion 
(AA) collisions vs. cm. energy per nucleon, ^/stm, measured by PHOBOS |2] (■), NA49 1 18] (★), and 
E895 iH (B) (see also ["^l); pp colHsions, by UA5 (A for non-single diffractive, T for inelastic events) at 
: 546 GeV |20| and ^/s^ = 200 and 900 GeV |21]; pp collisions (at lower ^/s^) from CERN-ISR 



(•) II22II and bubble chamber experiments 11231 12411 (SI) (the latter compiled and analysed in 12511 '). The 



pile 

inelastic UA5 data at y^s^ = 200 GeV is due to the extrapolation in f2\. The open symbols show the 



e+e^ measurements: the high-energy LEP mean multiplicities (o) averaged here from the data at LEP1.5 
= 130 GeV i2^|2j and LEP2 ^i^^ = 200 GeV |27, 28], and the lower-energy data by DELPHI 129] 
(□), TASSO (A), AMY (0), JADE (+), LENA (*), and MARKl (*) experiments. (See refs. in ] 30,,3liQM 
for e+e^ and pp/pp data). The solid line shows the calculations from Eq. (|3 based on our approach and 
using the corresponding fits (see text). The dashed and dotted lines show the fit |20] to the pp/pp data and 
the 3NL0 perturbative QCD 133] ALEPH fit l2Ji to e+e" data. The arrows show the LHC expectations. 

This nearness decreases the already small deficit in the e+e^ data as the energy 
increases. The deviation can be attributed to the inelasticity factor, or leading particle 
effect 1I14I1 in pp/pp collisions, which is known to decrease with the cm. energy. Then, 
at lower y/s^, some fraction of the energy of spectators contributes more into the 
formation of the initial state as the spectators pass by. This leads to the excess of the 
mean multiplicity in pp/pp data compared to the e+e data as it is seen in Fig. [T] 
Comparing further the average multiplicities from pp/pp collisions to those from AA 
ones, one finds that the data points are amazingly close to each other when the AA data 



are confronted the hadronic data at ^/s^ = 3 ^^'nn- The inclusion of the tripling energy 
factor indeed allows to describe such a fundamental variable as the mean multiplicity 
simultaneously in e+e^, pp/pp and central AA collisions for all energies. This shows that 
the multiparticle production process in headon AA collisions is derived by the energy 
deposited in the Lorentz-contracted volume by a single pair of effectively structureless 
nucleons similar to that in e+e annihilation and of quark-pair interactions in pp/pp 
collisions. Note that an examination of Fig. [T] reveals that not a factor 1/2 is needed 
to rescale ^/s^ to match the AA or e+e^ data as earlier was assumed for the mean 
multiplicity while recognised to unreasonably shift the e+e^ data on the pseudorapidity 
density at midrapidity when compared to the AA measurements [2]. This discrepancy 
finds its explanation in our consideration, within which the data on both the mean 
multiplicity and the midrapidity density (vide infra) are self-consistently matched for 
different reactions. Let us gain recall a factor 1/3 obtained earlier in lll6l] for ^/s^ for 
the pp mean multiplicity data relative to those from e+e data, similar to our finding. 

Fig.dshows that the mean multiplicities in different reactions are close starting from 
the SPS y^^NN , and become particularly close at ^/^nn ^ 50 GeV. However, at lower 
energies, the AA data are slightly below the e+e^ and hadronic data and the nuclear 
data increase faster with energy than the pp and e+e^ data do. On the other hand, 
as the cm. energy increases above a few tens GeV, the AA data start to overshoot 
the e+e^ data and reach the mean multiplicity values from pp interactions. From this 
one concludes on two different energy regions of the multiparticle production in AA 
reactions. The observations made can be understood in terms of the overlap zone and 
energy deposition by participants [10]. Due to this, one would expect the differences to 
be more pronounced in midrapidity densities as discussed below. 

3. In Fig.m we compare the pseudorapidity densities per participant pair at midrapid- 
ity as a function of ^/^nn from headon AA collisions at RHIC, CERN SPS and AGS to 
those of pp/pp data from CERN and Fermilab plotted vs. y^s^/3. Again one can see that 
up to the existing ^/s^ the data from hadronic and nuclear experiments are close to each 
other being consistent with our interpretation. The measurements from the two types of 
collisions coincide at 8 < ^/s^ < 20 GeV and are of the magnitude of the spread of 
AA data points at 200 GeV. However, above and below the 8-20 GeV region, there are 
visible differences in the midrapidity 7] -density values from AA vs. pp data. These in- 
dicate that, in contrast to the mean multiplicity which is a more global observable, the 
midrapidity density depends on some additional factor. As the densities are measured in 
the very central T] -region, where the participants longitudinal velocities are zeroed, it is 
natural to assume that this factor is related to the size of the Lorentz-contracted volume 
of the initial thermalized system determined by participants. 

To take into account the corresponding correction, let us consider our picture in the 
framework of the Landau model which is close to our description. Then, one finds for 
the ratio of the normalised charged particle rapidity density p{y) = {2/Npnn)dNch/dy at 
the midrapidity value = in AA reaction, Pnn» to the density Ppp in pp/pp interaction, 

Pnn(0)/Ppp(0) = 2 A^eh(^p/^NN) / (A^part A^^h") " (1) 

Here, A^^ch (^ch^ is the multiplicity in AA (pp/pp) collision, L = \n[y/s/(2m)], and m 
is the participant mass, e.g. the proton mass nip in AA reaction. According to our 
interpretation, we compare in the ratio ([T]) Pnn(O) to ppp(O) at ^/s^m = v^^/3 and 
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FIGURE 2. Pseudorapidity density p (0) of charged particles per participant pair (A^part / 2) at midrapid- 
ity as a function of the cm. energy of colHsion. The open and combined symbols show the pseudorapidity 
density values vs. cm. energy per nucleon, -^snn, measured in the headon AA collisions by BRAHMS jal 
(o), PHENIX |7] (A), PHOBOS |2] (□), and STAR |341 (★), and the density values recalculated ff\ from 
the measurements taken by CERES/NA45 |35] (+), NA49 |36| (0), NA50 |37] (o) WA98 |38] (*), E802, 
and E917 |39] (*). The nuclear data at ^^s^ around 20 GeV and the RHIC data at y/I^ = 130 GeV 
and 200 GeV are given spread horizontally for clarity. The solid symbols show the pseudorapidity density 
values vs. cm. energy as measured in non-single diffractive pp collisions by UAl |401 (■), UA5 

I?, 201 (A), CDF isj] (▼), and from inelastic pp data from ISR |22] (★), and bubble chamber |24, 41] (•) 
experiments (the latter as recalculated in |4]). The solid line connects the predictions from Eq. (|2j. The 
dashed line gives the fit to the calculations using the 2nd order log-polynomial fit function analogous to 
that used 1 5 ] in pp data. The fit function from 1 5] is shown by the dashed-dotted line. The dotted line shows 
the linear log approximation of UA5 to inelastic events |4]. The arrows show the LHC expectations. Note 
that e+e^ data at y^i^ = 14 GeV to 200 GeV (not shown) follows the heavy-ion data 01 • 

consider a constituent quark of mass ^nip as a participant in pp/pp collisions and a proton 
as an effectively structureless participant in headon AA collisions. Then, Eq. ([T]) reads: 

(2) 



Pnn(O) = 2A^chPpp(0) Jl -41n3/ln(4m2/5NN) /(A^partA^Ph^) • 



Using the fact that the transformation factor from ytorj does not influence the above 
ratio and substituting the multiplicity values from Fig.[T]and of ppp(O) from Fig.|2linto 
Eq. ©, one obtains the values of Pnn(O), displayed in Fig.Elby solid line. One can see 
that the correction made provides good agreement between the calculated Pnn(O) values 



and the data. Eq. © shows the importance of the correction for the participant type to be 
introduced as argued above. One can see that our calculations account also for different 
types of rise of AA data below and above SPS region. Note that the same two regions 
recently have been indicated by PHENIX CD from the ratio of the midrapidity transverse 
energy density to the pseudorapidity density. From these findings, one can expect the 
midrapidity transverse energy densities in pp/pp and headon AA collisions to be similar 
due to the description proposed here. Also, the SPS transition region properties discussed 
byNA49 0], can be treated without any additional assumptions. 

4. To estimate Pnn(O) for aV^nn > 200 GeV, we extrapolated the values of Eq. ^ 
utilizing the function found fy\ to fit well the pp data. The predictions for Pnn(O) and 
the fit for pp data are shown in Fig.|3by dashed and dashed-dotted lines, respectively. 

The obtained Pnn(O) show faster rise with ^Js^^ than Ppp(O). Our calculations, 
sharing the behaviour at SPS-RHIC enermes with that up to the LHC ones, give 
Pnn(O) ~ 7.7 for LHC. From the CDF fit [5] and assuming it covers LHC energies, 
one finds ppp(O) ^ 6. 1. Our Pnn(O) value for LHC is consistent with that of ^ 6. 1 given 
in the PHENIX extrapolation [7] within 1-2 particle error acceptable in the calculations 
we made. Our result is in a good agreement with the best ATLAS Monte Carlo tune 
1E2I1 . Noticing that ^/s^ is near to ^/5^/3 at LHC, the close values of Pnn(O) and 
Ppp(O), predicted for LHC by us and estimated independently in [i2i,i42], demonstrates 
experimentally grounded description and predictive ability of our interpretation. 

Solving Eq. (0) for A/ch/( 0.5 A/part) we predict the A A mean multiplicity energy de- 
pendence at ^/s^ > 200 GeV. In this calculations, we use the fits of Ppp(O) fs"] and A'^^j^ 
0] and our approximation for Pnn(O), all shown in Figs. [T] and |2l From the resulted 
curve for A^ch/ (O.SA^pan) given in Fig.[Tl one finds that the value obtained for LHC is just 
about 10% above the N^^iV^) fit [20] prediction for LHC and about 3.3 times larger 



the AA RHIC data at ^^nn = 200 GeV. Again, this number is comparable with the es- 
timate made by |7] and points out to no evidence for change to another regime as the 
^5nn increases by about two magnitudes from the top SPS energy. Nevertheless, one 
can see that the data obtained at the highest RHIC energy give a hint to some border-like 
behaviour of the mean multiplicity where the pp/pp data saturate the nuclear data, and 
another transition energy region is possible to be found (as at low energies). This makes 
AA experiments at ^^nn > 200 GeV of particular interest. 
5. At the end, let us dwell on the following. 

From our description, the mean multiplicity in nMcfeon-nucleus collisions is predicted 
to be of the same values as that in pp/pp data, and, moreover, almost no centrality 
dependence is expected for such type of interactions These predictions are well 
confirmed by various data from hadron-nucleus collisions at ^/^nn ~ 10-20 GeV to 
recent RHIC dAu data at 200 GeV 10] • The same seems to be correct also for the 
pseudorapidity density at midrapidity, which is already supported to be a trend ||2(]. These 
findings remind about similar conclusions made about two decades ago fl2!]. 

The recent observation [8, 9] made at RHIC for multihadron data from CuCu colli- 
sions to not change compared to same cm. energy AuAu data when scaled for the same 
participant numbers is also understood due to our description as already mentioned. In- 
deed, for the same number of participants, no difference in the bulk variables is expected 
as one moves from one type of (identical) colliding nuclei to another one at the same 



cm. energy as soon as the same energy is deposited into the thermalization zone. Note 
that the proper definition of participants and, thus, of the energy available for particle 
production^ we discuss here, allows scaling within the constituent quark picture to be 
applied @,Ii3>S3] to model the multihadron data at RHIC for different observables. 
One of us (EKGS) is grateful to Organizers for invitation and partial financial support. 
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